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ABSTRACT 



Various theoretical and observational results have been reported regarding 
the presence/absence of net electric currents in the sunspots. The limited spa- 
tial resolution of the earlier observations perhaps obscured the conclusions. We 
have analyzed 12 sunspots observed from Hinode (SOT/SP) to clarify the issue. 
The azimuthal and radial components of magnetic fields and currents have been 
derived. The azimuthal component of the magnetic field of sunspots is found 
to vary in sign with azimuth. The radial component of the field also varies in 
magnitude with azimuth. While the latter pattern is a confirmation of the in- 
terlocking combed structure of penumbral filaments, the former pattern shows 
that the penumbra is made up of a "curly interlocking combed" magnetic field. 
The azimuthally averaged azimuthal component is seen to decline much faster 
than 1/vj m the penumbra, after an initial increase in the umbra, for all the 
spots studied. This confirms the confinemen t of magnetic fields and absence of 
a net current for sunspots as postulated by iParkerl (119961 ). The existence of a 
global twist for a sunspot even in the absence of a net current is consistent with 
a fibril-bundle structure of the sunspot magnetic fields. 



Subject headings: Sun: magnetic fields. Sun: photosphere. Sun: sunspots 



1. Introduction 



Sunspots have shown evidence for twist even from the time of iHald (119251 . 119271 ) who 
postulated the hemispheric rule for the chirality of chromospheric whirls. This was later 
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confirmed with a larger data set by iRichardsoru (1194 ll ). Evidence for pfiotosplieric cliiral- 
ity could be seen in early continuum images of sunspots, obtained with exceptional image 
quality. Later, photospheric vector magnetogram s showed global twis t inferred from the 
non- v a nishing average s of the force- free parameter (jPevtsov et al.l (Il994j ) ; iHagino &: Sakurai 
( I2OO4J ): iNandyl (120061 ) and references the rein). The rion- force -free nature of photospheric 
magnetic field in the sunspots, prompted iTiwari et al.l (l2009al ) to propose the signed shear 
angle (SSA) as a more robust measure of the global twist of the sunspot magnetic field. 

Although, the sign of SSA matches well with the sign of the global alpha parameter, 
the magnitudes are not so well correlated. The physical significance of a globally averaged 
a parameter rests heavily on the existence of a net current in the photospheric sunspot 
magnetic field. One way of arriving at a global a is by taking the ratio of total vertical 
current to the total flux ('in tegral method). This va lue was found to agree with the values 
obtained by other methods (IHagino fc Sakuraill2004l ). 



For a monolithic sunspot magnetic field, the global twist and net current is expected to 
be well correlated by Ampere's Law. However, the existence of a net current is ruled out 
theoretically for fibril bundles as well as for monolithic fields with azir nuthal field d ecreasing 
faster than l/cc, where w is the radial distance from the spot center (jParkerlll996l ). Several 
attem pts to resolve this problem using vector magnetograms ha ve not been very conclusive 
so far Jwilkinson et al.lll992l : iLeka et al.lE996l : Iwheatlandlboooh . 



A resolution of this problem can be used to disentangle the relation between global 
twist and the global a parameter. Also, the resolution is needed to evaluate the so called 
hemispheric hel icity rule seen in the global a parameter calculated from photosp heric vector 
magnetograms (jPevtsov et al.lll994l . Il995l : IHagino fc Sakuraill2004l : iNandyl 120061 ). The avail- 
ability of high resolution vector magnetograms from Hinode (SOT/SP), gives us the best 
opportunity so far to address this problem. The effec t of polarimetric noi se is expected to be 
negligible in the estimation of magnetic parameters (jTiwari et al. 2009bl ) from these data. 



In this Letter we obtain an exp ression for the net current using a generalization of the 
expression obtained by iParkerl (119961 ). We then proceed to measure this current from several 
vector magnetograms of nearly circular sunspots. We finally discuss the results and present 
our conclusions. 



2. Expression For Net Current 



Following iParkerl (119961 ) , we consider a long straight fiux bundle surrounded by a region 
of field free plasma. We use the words "field free" in the empirical sense that there is no large 
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scale coherent and unipolar magnetic field surrounding the flux b undle. Also, w e include 
the case where the bundle can be replaced by a monolithic field. iParkerl (119961 ) assumed 
azimuthal symmetry as well as zero radial component B^., of the magnetic field. For realistic 
sunspot fields, we need to relax both these assumptions. 

The vertical component of the electric current density consists of two terms, viz. — —^Mf- 



term as the "twist current density" 
center is then given by 

iJw) 



Jt- 



fJ-or dtp 

second 

The net current within a distance w from the 



and We will call the first term as the "pleat current density", j„ and the second 

^or or f J ^ jp 
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The ip integral over jp vanishes, while the second term yields 



w 



(2) 



which gives the net current within a circular region of radius w. 



3. The Data Sets and Analysis 



We have analyzed the vector magnetogra r ns obtained frorn Sola r Optical Telescop e /Spec tro- 
polarimeter (SOT/SP : iTsuneta et al. IShimizu et al.l (120081 ): ISuematsu et al.l (120081 ) 
Ichimoto et al.l (120081 )) onboard Hinode (IKosugi et al.l 120071 ). The calibration of data sets 
have been performed using the standard "SP_PREP" routine developed by B. Lites and avail- 
able in the Solar- Soft package. The prepared polarization spectra have been inverted to ob- 
tain v ector magnetic field components using an Unno-Rachkowsky ( Unno 19561: Rachkowskv 



1982 



1967 ) inversion under the assum ption of Milne-Eddington (ME) atmosphere (iLandolfi &: Landi Degl'Innoce: 



Skumanich &: Liteslll987l ). We have used the inversion code "STOKESFIT" which has 



been kindly made available by T. R. Metcalf as a part of the Solar-Soft package. We have 
used the newest version of this code which returns true field strengths along with the filling 
factor. The azimuth determination has inherent 180° ambiguity due to insensitivity of Zee- 
man effect to orientation of the transverse fields. Numerous techniques have been developed 
and applied to resolve this problem, but not even one guarantees a complete resoluti on. The 
180° azimuthal ambiguity in our data sets are rem oved by using acute angle method (IHarvey 
19691 : ISakurai et al.l 1 19851 : ICuperman et al.lll992l ). 



In order to minimize the noise, pixels with transverse [Bt) and longitudinal magnetic 
field (Bz) greater than a certain level are only analyzed. A quiet Sun region is selected for 
each sunspot and Icr standard deviation in the three vector field components B^, By and 
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Bz are evaluated separately. The resultant standard deviations of and By is then taken 
as the la noise level for transverse field components. Only those pixels where longitudinal 
and transverse fields are simultaneously greater than twice the above mentioned noise levels 
are analyzed. The data sets with their observation details are given in Table 1. We have 
treated each polarity as an individual sunspot whenever both the polarities are observed 
and compact enough to be studied. We have studied only those spots where the polarity 
inversion lines are well separated from the edge of the sunspot. 

The results of the inversions yield the 3 magnetic parameters, viz. the field strength 
B, the inclination to the line of sight 7, and the azimuth 0. These parameters are used to 
obtain the 3 components of magnetic field in Cartesian geometry as 

5^ = 5 cos 7 (3) 
By = 5 sin 7 sin (4) 
Bx = B sin 7 cos (5) 



This vector field is transformed to heliographic coordinates ( IVenkatakrishnan fc Garylll989l ) 
for the spots observed at viewing angle more than 10°. The transverse vector is then ex- 
pressed in cylindrical geometry as 

Br = -{xBx + yBy) (6) 
r 

B^ = ^{-yBx + xBy) (7) 

The azimuthal field B^ is then used in equation (2) for obtaining the value for the total 
vertical current within a radius w. 

We have computed "twist angle" for all the sunspots using B^ and B^. as shown in Table 
1. The error in "twist" measurement is simply the error in azimuth measurement. Using the 
weak field approximation, we can find the azimuth -0 from tan 2-0 = U /Q. From this we can 
estimate the error in if) as equal to the percentage error in linear polarization measurements. 
Thus, a 1% error in polarimetry means that the error in if) equals 0.01 radians or 0.57 degrees. 
We have performed Monte Carlo simulations of the effect of noise on the inversions which we 
plan to present in a more detailed paper. We have verified that the error in ip is consistent 
with the value estimated from the weak field approximation. 



We can see in the Table 1 that the twist angles for regular sunspots match well with 
the global SSA as expected, whereas they do not match for irregular sunspots. 
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The Results 



Figure 1 shows an example of the maps of twist current, pleat current, and Br 
for a sunspot NOAA AR 10933 which is nearly circular. Figure 2(a) shows plots of B^ 
and Br along with the different concentric circles around spot center. In Figure 2(b), the 
spatial variation of both B^ and Br are clearly seen. This variation is corresponding to a 
typical circle selected in the penumbra. The Br variation in the penumbra is a ma nifestation 
of the interlocking combed structure (jlchimoto et al.l 120071 : iTiwari et al.l l2009al ). The B^ 
variation in the penumbra shows that not only is there an interlocking combed structure, 
but these structures are curled as well. In other words, we may describe the penumbral field 
as possessing a "curly interlocking combed" structure. This fe ature of the deviation of the 
vector field azimuths from a radial direction was also seen by iMathew et al.l (120031 ) in the 
magnetic field of a sunspot belonging to NOAA AR 8706, using the infra-red Fel line pair 
at 1.56 micron. 

The azimuthal averages < B^ > and < Br > were obtained at different values of 
vj. Figure 2(c) shows the plots of < B^ > and < 5^ > as a function of vj. The circles 
corresponding to the selected radii are shown in the upper panel of the same figure. The 
azimuth- averaged < Br > drops rapidly to a very low value at the edge of the sunspot. This 
is a clear evidence for the existence of a canopy where the field lines lift up above the line 
forming region. Figure 3 shows the plot of log < B^ > as a function of log w. The slope 6 
of the declining portion of this plot is 9.584, which shows that field varies faster than 1/w. 
This can be construed as evidence for the neutralization of the net current. The 6 for other 
sunspots have also been computed and are given in Table 1. 

The map of vertical current density for the same sunspot is shown with intensity scale 
in the left panel of Figure 4. The values are expressed in Giga Amperes per square meter 
(GA/m^). We can see that the dis tribution of u is dora inated by high amplitude fluctuations 
on small scale as also reported in lTiwari et al.l (l2009al ). It is therefore difficult to make out 
any systematic behaviour of the sign of as a function of w. 

The right panel of the same figure shows the total current within a radius w as a function 
of zu. As expected from the trend in Figure 3, the total current shows evidence for a rapid 
decline after reaching a maximum. Similar trends were seen in other sunspots. We have also 
plotted in right panel of Figure 4, the net current as calculated by the derivative method 
(viz. summation of current densities calculated as the local curl of B) . We do see a trend of 
neutralization, although the effect is less pronounced because of the larger noise present in 
the derivative method. We can also infer from the right panel of Figure 4 that the increments 
of net vertical current flowing through annular portions of the sunspot do show a reversal in 
sign. 
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Table 1 shows the summary of results for all the sunspots analyzed. Along with the 
power law index 6 of decrease, we have also shown the average deviation of the azimuth 
from the radial direction ("twist angle"), as well as the SSA. The average deviation of the 
azimuth is well correlated with the SSA for nearly circular sunspots, but is not correlated 
with SSA for more irregularly shaped sunspots. Thus, SSA is a more general measure of the 
global twist of sunspots, irrespective of their shape. 



5. Discussion and Conclusions 



It is well known for astrophysical plasmas, that the plasma dist orts the mag netic field 
and the curl of this distorted field produces a current by Ampere's law (jParkerlll979l ). Parker's 
(1996) expectation of net zero current in a sunspot was chiefly motivated by the concept of 
a fibril structure for the sunspot field. However, he also did not rule out the possibility of 
vanishing net current for a monolithic field where the azimuthal component of the vector field 
in a cylindrical geometry declines faster than While it is difficult to detect fibrils using 
the Zeeman effect notwithstanding the superior resolution of SOT on Hinode, the stability 
and accuracy of the measurements have allowed us to detect the faster than 1/m decline of 
the azimuthal component of the magnetic field, which in turn can be construed as evidence 
for the confinement of the sunspot field by the external plasma. The resulting pattern of 
curl B appears as a drop in net current at the sunspot boundary. 

If this lack of net current turns out to be a general feature of sunspot magnetic fields 
in the photosphere, then measurement of helicity from a global average of the force-free 
parameter becomes suspect. On the other hand, sunspots are evidently twisted at pho- 
tospheric levels, as seen from the non-vanishing average twist angle as well as the SSA 
(Table 1). Although the existence o f a global t A vist in the absence of a net current is pos- 
sible for a monolith ic sunspot field (iBatyl 120001 : lArchontis et al.l l2004l : iFan &: GibsonI |2004J : 



Aulanier et al.ll2005l ). a fibril model of the sunspot field can accommodate a global twist even 
without a net current (jParkerlll996l ). 



The spatial pattern of current density in a sunspot (e.g., left panel of Figure 4) is really 
a manifestation of the deformation of the magnetic field (V x B) by the forces applied by the 
plasma. The Lorentz force exerted by the field on the plasma produces an equal and opposite 
force by the plasma, thereby confining the field. Thus our analysis actually shows the pattern 
of the forces exerted by the plasma on the field. The sharp decline of the azimuthal field 
with radial distance thus shows the confinement of the sunspot magnetic field by the radial 
gradient of the plasma pressure. 
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Theoretical understanding of the penumbral fine structure has improved considerably 
in recent times (IThomas et al.ll2002l : IWeiss et aLll2004l ). The onset of a con yective ins t abihty 
for magnetic field inchna tion exceeding a critical valu e was proposed by iTildeslevI (120031 ) 
and lHurlburt et al.l (120021 ) . A bifurcation in the onset (IRucklidge et al.lll995l ) could explain 
other features like hysteresis in the appearance of penumbra as a f unction of sunspot size . 
Numerical sirnulation of magneto-convection also steadily improved (IHeinemann et al.ll2007l : 



Rempel et al 



( iRempel et al. 



2009a[) . culminating in very realistic production of penumbral field structure 



2009bl ). It is possible, owing to the random and stochastic nature of convective 



structures, that no net twist in the simulated spot field would be produced by convection for 
negligible Coriolis force. If so, it would be very interesting to simulate magneto-convection 
in a twisted sunspot field. In this case, would the resulting fine structure mimic the observed 
"curly interlocking combed" structure of the penumbral magnetic field? If not, we must look 
elsewhere for explaining the "curly interlocking combed" structure. A twisted fibril bundle 
would then be a sol ution. Recent examples of filamentary penurnbral structures based on 
such cluster models ( ISolanki fc Montavonlll993l : ISpruit fc Scharmeiil2006l : IScharmer fc Spruit 
20061 ) have also been proposed. 



Parker! (119961 ) also mentions the possibility of net currents in the corona, continuing 



down to the height where the first cleaving takes place. It would therefore be imperative 
to look for net currents at highe r reaches of th e solar atmosphere. This is very important 



because several theor ies of flares (jMelrosdll995l ) and CME triggers (IForbes fc Isenberg||l991 



Kliem fc Torokl 2006) rely heavily on the existence of net currents in the corona above the 
sunspots. 

Future large ground based telescopes equipped with adaptive optics and multi spectral 
line capabilities would go a long way in addressing these issues. In the meantime, direct 
measurement of the global twist of sunspots using parameters like the SSA should serve as 
proxies for estimating the net currents of active regions in the corona. The SSA will also be 
a better parameter to base a fresh look at the hemispheric rule in photospheric chirality. 
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Tablc 1: List of the active regions studied. The power index 6: the slope of decrease of B,!, 
value, the twist angle, the signed shear angle (SSA) and other details of the sunspots are 



given: 



AR No. 


Date of 


Slope 


Shear Angle 


Twist Angle 


Position 


Hemispheric 


(NOAA) 


Observation 


S 


(SSA: deg) 


(tan-HB^/Br): deg) 




Helicity Rule 


10969 


29 Aug 2007 


7.514 


-4.488 


-4.009 


S05W33(t) 


No 


10966 


07 Aug 2007 


4.349 


-5.120 


-7.028 


S06E20(t) 


No 


10963(-) 


12 Jul 2007 


4.366 


-5.123 


41.8731' 


S06E14(t) 


No 


10963(+) 


12 Jul 2007 


4.210 


-4.495 


-5.112 


S06E14(t) 


No 


10961 


02 Jul 2007 


4.976 


-4.973 


29.451" 


S10W16(t) 


No 


10960 


07 Jun 2007 


3.267 


3.182 


-24.012" 


S07W03 


Yes 


10953 


29 Apr 2007 


8.249 


-3.382 


7.200" 


S10E22(t) 


No 


10944 


03 Mar 2007 


2.407 


-4.635 


-5.130 


S05W30(t) 


No 


10940 


01 Feb 2007 


2.281 


-4.726 


-7.950 


S04W05 


No 


10933 


05 Jan 2007 


9.584 


-2.283 


-2.689 


S04W01 


No 


10926 


03 Dec 2006 


2.750 


-1.538 


6.001" 


S09W32(t) 


No 


10923 


10 Nov 2006 


3.175 


0.785 


-9.010" 


S05W30(t) 


Yes 



Note. — 1. (t) : transformed 
— 2. ti •■ twist angle for irregular sunspots does not fit to a cylindrical assumption and therefore gives incorrect 
values. 
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Fig. 1. — Examples of the two components of the vertical electric current density namely 
the "twist" and the "pleat" current densities {jt and jp) observed in NOAA AR 10933 are 
shown in the upper panel. The lower panel shows the azimuthal and radial component of 
the magnetic field {B^ and B^). 
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Fig. 2. — (a): The same image as the lower panel of Figure 1 but with concentric circles over 
plotted on them, (b): Plots of and Br along the periphery of a typical circle (45*'^ pixel 
away from center) selected in the sunspot. (c): The mean B^ and mean Bj. with ^ and | 
of their variations respectively along radial direction with different zu have been plotted. 
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Fig. 3. — Variation of log < > with log to. The slope of the declining portion of plot has 
been calculated as a simple power law index 6 for each sunspot and has been given in Table 
1. 
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Fig. 4. — Left panel: The map of vertical current density is shown with intensity scale. 
The values are expressed in Giga Amperes per square meter (GA/m^). Right panel: The 
net current variation with increasing area has been shown. The solid line shows the results 
of the calculations from the Equation 2. Also shown, by a dashed line, is the results from 
the derivative method. We can see the net current reduces very fast after a peak and 
almost vanishes for complete sunspot. On the other hand the net current computed from 
the derivative method shows a shallow behaviour. 



